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I. INTRODUCTION 



Spin gap systems, where the spin gap is of quantum 
mechanical origin, are of interest to both theoretical and 
experimental investigations. These systems have spin- 
disordered ground states which can be described as quan- 
tum spin liquids^. Spin gap is known to open, for exam- 
ple, in 5 = 1/2 spin systems that form lattices of coupled 
dimers^, giving rise to many interesting phenomena. The 
Shastry-Sutherland modcl^ provides a unique example 
of such two-dimensional frustrated networks of S" = 1/2 
dimers. The model includes nearest (J) and next nearest 
neighbor (J') antiferromagnetic interactions as shown in 
Fig. [TJ with the Hamiltonian 
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FIG. 1. (color onUne) The Shastry-Sutherland lattice in the 
SrCu2 (303)2. The orthogonal dimers (thick lines) are de- 
noted by letters A and B, while the Cu ions on the dimer are 
enumerated as 1 and 2. The inter-dimer bonds (with Heisen- 
berg exchange J') are shown by a thin line. We also show 
the Si and C2v symmetry point groups of the buckled CuBOs 
layer (the open and closed circles indicate the Cu ions that 
are below or above the layer, respectively). 



In the case of J' = the model is reduced to a lattice 
of independent dimers, where in the ground state the 
S = 1/2 spins of each dimer form a singlet, and the 
ground state wave function is just the product of these 
independent dimer-singlcts. According to Shastry and 
Sutherland^, the singlet dimer product state is an exact 
eigenstate of the Hamiltonian even for finite values of the 
J' due to the particular geometry of the lattice. 

An experimental realization of the Shastry-Sutherland 
model is the quasi two-dimensional antiferromagnetic 
compound^ SrCu2 (603)2- This material has a tetrag- 
onal unit cell and is characterized by the alternating lay- 
ers of CuBOa molecules and Sr^''' ions; in the former, the 
Cu^+ ions occupying crystallographically equivalent sites 
carry spin S = 1/2 degrees of freedom and form a lattice 
of orthogonal dimers connected by the triangular shaped 
BO3 molecules^i^. A schematic figure of CuBOs layer is 
shown in Fig. [1] 

Magnetic susceptibility measurements of this mate- 
rial revealed a peak at around 20 K and a sharp drop 
to zero at decreasing temperatures^. Fitting the ex- 
ponential curve that is characteristic for the spin gap 
systems, Kagcyama et al. estimated in Ref. [1] the 
gap to be A w 19 K, while from the NMR relaxation 
rate they obtained the gap about 30 K. Magnetization 
measurements^ clarified the presence of a gapped spin- 
singlet ground state and a continuos transition to the 
gapless magnetic state at 20 T which corresponds to a 
gap of 30 K in a good agreement with the relaxation 
rate measurements. While early magnetization measure- 
ments in high fields revealed plateaus only at 1/4 and 1/8 
of the saturated magnetization^, refined measurements 
have suggested^ more plateaus at 1/3 and other values 
of magnetization. 

Miyahara and Ueda^ pointed out that SrCu2 (603)2 
can be satisfyingly described by the Shastry-Sutherland 
model. They estimated the critical point where the sin- 
glet dimer ground state goes to the Neel-state to be 
{J'/J)c ~ 0.7 performing variational calculations and 
exact numerical diagonalization (on the basis of scries 
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expansion^ and exact diagonalizationiS, it is now be- 
lieved that the above two states are mediated by a new 
plaquette-singlet phase and that a transition from the 
dimer phase to the plaquette-singlet occurs at {J'/J)c = 
0.68). Using the experimental findings of Ref. [Jl they 
estimated the nearest neighbor coupling constant to be 
J = 100 K and the next nearest neighbor coupling J' = 
68 K. This yields J'/ J = 0.68 indicating SrCua (803)2 to 
be close to the transition point ( J'/ J)c = 0.7. Later this 
estimate has been updated to J = 7.3 meV with J' / J = 
0.635 (Ref. [ill) or J = 6.16 meV with J'/J = 0.603 
(Ref. M)- 

Furthermore, Miyahara and Ueda carried out pertur- 
bation theory in the dimer state up to the fourth order 
in J'/J and found that the triplet excitations are local- 
ized. The hopping of triplets is only possible through 
closed paths of dimer bonds, thus only from the sixth 
order in perturbation. This property of the triplet exci- 
tations is related to the formation of plateaus. At certain 
values of the magnetization the excitations localize into 
a supcrlattice structure to minimize the energy^. Mo- 
moi and Totsuka— have explained the appearance of 
plateau states through the scenario of metal to Mott- 
insulator transition where the triplet excitation were 
treated as bosons interacting via various repulsive in- 
teractions arising from higher-order perturbation in J'. 
At dominating repulsive interaction, the triplet exci- 
tations localize and crystallizes in commensurate pat- 
terns developing the plateau states. In fact, NMR spec- 
troscopy by Kodama et al. exhibitedi^ directly the su- 
perlattice structure at m/msat = 1/8. Recently new 
magnetization plateaus have been foundiS by nonpertur- 
bative Contractor-Renormalization (CORE) method at 
1/9, 1/6 and 2/9 of the saturation, while the analysisii 
using the pcrturbative continuous unitary transformation 
(PCUT) has predicted, on top of the above ones, one 
more plateau at 2/15. It has also been argued that that 
the inclusion of the spin-lattice effects determines the 
spin structure in the plateaus.— 

In the past few years various experiments have been 
carried out to examine the spin excitations. While 
the original Shastry- Sutherland model is isotropic in 
spin space, its experimental realization SrCu2 (603)2 ex- 
hibits anisotropic behavior; inelastic neutron scattering 
mcasurcmentSii^, electron spin resonance^S, and Raman 
scattering^! indicated a splitting of the triplet excita- 
tions at the F-point, which was explained to be caused 
by the effect of inter-dimer Dzyaloshinskii-Moriya (DM) 
vector directed perpendicular to the copper plane^^. 
Later another splitting was found^^ at the q = (tt, 0) 
point indicating the relevance of in-plane components 
of the DM interaction. The ESR study of Nojiri et 
al.— shows an anti-level crossing at the critical mag- 
netic field where the lowest-lying triplet excitation would 
cross the singlet level, which is consistent with the per- 
sistent spin gap found in the specific-heat^i- and the 
NMR measurements^^. These splittings and the anti- 
level crossing mean that states corresponding to different 



magnetization (singlets and triplets) are mixing and 
is no longer a good quantum number. This mixing be- 
tween the singlet and triplet states of a dimer can be ex- 
plained by an intra-dimer anisotropy, e.g. an intra-dimer 
DM vector. For more details on the Shastry-Sutherland 
model and SrCu2 (603)2, we refer the readers to the re- 
view articles Ref. (theory) and Ref. [l^ (experiments) . 

Magnetization process in such dimer systems as 
TlCuCla is fairly well understood; onset of magnetiza- 
tion is triggered by Bose-Einstein condensation of gapped 
triplons and the magnetic phase above the critical field is 
characterized by broken XY-symmetry perpendicular to 
the applied field (see, e.g., Ref. [28l for a review). Dynam- 
ics at high-fields is also well described within the above 
scenario^^. On the other hand, the existence of DM in- 
teractions is known to substantially modify the above 
picture and even new phases may appear in the presence 
of DM interactions. Moreover, the small kinetic energy 
and relatively large (effective) interactions among dimcrs 
lead to various magnetic superstructuresi^i^, which we 
cannot simply neglect in considering the dynamics at 
high fields. Nevertheless, the global structure of the 
phase diagram and the dynamics in the presence of DM 
interactions and magnetic superstructures is only par- 
tially understood. The aim of this paper is to present a 
simple theoretical framework to investigate the ground- 
state phases and the magnetic excitations over them with 
the extension to the cases with superstructures in mind. 
Specifically, by using the bond-wave approximation, we 
examine the excitation spectrum of SrCu2 (603)2 at zero- 
and low magnetic fields below the plateaus and com- 
pare the results with the neutron-scattering- and the ESR 
data. 

The paper is structured as follows: in Sec. II we review 
the symmetry group of SrCu2 (603)2 and determine the 
allowed anisotropic terms in the Hamiltonian. In Sec. 
Ill we describe the variational approach we use to get 
the ground state and the way we construct the excitation 
spectrum using the bond operators. In Sec. IV, we show 
the variational phase diagram and the spectrum for the 
zero field case. In Sec. V we map out the phase diagram 
in the presence of a field perpendicular to the basal plane. 
In Sections VI and VII we describe the ESR spectra in 
magnetic fields perpendicular and parallel to the plane, 
respectively. Last, we compare our results with neutron 
scattering experiments and the ESR spectra in Sec. VIII. 
We conclude with Sec. IX. 



II. SYMMETRY CONSIDERATIONS AND THE 
MODEL HAMILTONIAN 

At high temperatures the spacegroup of SrCu2 (603)2 
is a tetragonal lA/mcm (Refs. and [3l|). A struc- 
tural distortion in the CU6O3 layers below Tg = 395 K 
shifts the two types of orthogonal dimer planes along 
the z axis in opposite direction, lowering the symmetry 
of the SrCu2 (603)2 to /42m at low temperatures, also 
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FIG. 2. (color online) The components of the symmetry 
allowed intra-dimer (a) and inter-dimer (b) Dzyaloshinskii- 
Moriya vectors shown in a unit cell. The arrows running from 
i to j on the bonds indicate the ordering of the spin operators 
Si X Sj in the DM interaction term. 



with tetragonal symmetr}*^. Restricting ourselves to the 
symmetries of the CuBOs layer, above T, the two — to 
each other orthogonal — types of dimers lay in the same 
plane and the wallpaper group pAg consists of the point 
group = {E,C2{z),axz,o'yz} at the middle of dimer 
bond and of the point group C4/1 in the middle of four 
sites belonging to different dimers. Lowering the temper- 
ature, while the C2V remains a symmetry of the buckled 
CuBOs layers, the loss of the ah reflection plane below 
Tg lowers the Cih point group to 54 = {£', ^4, 6*2(2 ), 5*1} 
(see Fig. [T]).The unit cell in both cases consists of two 
orthogonal dimers — dimer A that is parallel to the x 
axis, and dimer B parallel to y axis. 

The symmetry of the lattice determines the possible 
terms in the Hamiltonian: the components of the 5-tensor 
anisotropy and of the exchange interactions, including 
the components of the DM interactions; see Table Hand, 
e.g., Ref. [25|. Here we shortly present the relevant terms. 
The sites Al, A2, Bl, and B2 correspond to the sites 1, 
2, 3, and 4 of Ref. [23, respectively. 

The 5-tensor at site Al takes the following form: 

f 9x -gs \ 
gAi = gy (2) 
\-9s g, J 

with gx — gy as required by the tetragonal symmetry. 
When the field is in the z direction, the Zeeman term 
reads 

Hh = - gzmHz (S'ai + 5'a2 + -Sm + 8^2) ^^^^ 
+ gsfJ-BHz (S'ai - 5'a2 + -Sgi - -532) ' 



while if the field is along the x axis, it reads 

= - gxV-BHx (5ai + 5*12) - 9yHx (<5|i + S'12) ^ 
+ gsMBi/x(5ii-5i2) . ^ ^ 

For convenience, we choose hz = gz^-sHz and introduce 
the scaled variable gs = gs/gz- 

Next, we consider the DM interactions Hdm = Hd + 
% D' , with the intra-dimer 

Hz5 = ;^D,,-(S, X S,) (4) 

NN 

and inter-dimer 

Hz,. = ^D',,.(S, xS,) (5) 

NNN 

contributions. The summation is over the nearest neigh- 
bor (NN) and next-nearest neighbor (NNN) sites. Once 
we specify the DM vectors on a bond, the DM interac- 
tions on the remaining bonds of the unit cell follow from 
the symmetry of the cell, as shown in Fig. [51 To specify 
the sign of the DM interaction unambiguously, we have 
also denoted by an arrow the order of the spins in the 
cross product: an arrow from site i to site j in the figure 
means that we need to take Si x Sj. 

At temperatures above the structural transitior*^ T > 
Tg = 395 K (high-temperature phase), the middle of a 
bond is an inversion center due to which there is no DM 
vector on the dimer bond and only the inter-dimer D' 
perpendicular to the CuBOa plane is allowed: 

D = , (6a) 
D' = (0,0,i:'l) (6b) 

(we denote by I?^ the z component). In the following, 
we will call this case the high-symmetry case. 

Below Ts (low temperature phase), however, this in- 
version symmetry is lost and the in-plane DM compo- 
nents are allowed as well. Correspondingly, D' becomes 
an arbitrary vector and the intra-dimer D is lying in the 
CuBOs-plane and perpendicular to the dimer, so that 

Da = (0,Z?,0), (7a) 
Db = (-Z?,0,0) , (7b) 

with the site ordering convention as shown in Fig. [21 We 
refer to this case as the low-symmetry case. 

A discussion of the different estimates for the strength 
of the DM interactions and the p-tensor anisotropics for 
SrCu2(B03)2 is presented in Sec. lVlIIl 

III. THE VARIATIONAL APPROACH AND 
THE BOND-WAVE THEORY 

A. Variational wave function 

The ground state of the pure Shastry-Suthcrland 
model (|T]) can be written as a product of singlets \s) over 
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TABLE I. Summary of symmetry analyses. If the g-tensor 
anisotropy is taken into account, the spin 0(2)-symmetry is 
lost. 

high-symmetry low-symmetry 

symmetry (unit cell) T>4h 2?2ii 

DM (intra-dimer) forbidden D || (a6) A D _L dimer 

DM (inter-dimer) D' || c arbitrary 

spin (h II z) 0(2)-sym — 

spin {h \\ x) — — 



the dimer bonds, 5' = ridimcrs I*)- presence of the 

DM interactions and finite magnetic fields, we need to 
extend this wave functions to a variational one. Namely, 
we allow for a linear combination of the singlet and triplet 
states on each dimer (we keep the dimer wave function 
entangled), while we retain the product form over the 
dimer bonds: 



I*) = 



n i^A> n i^'i 



A dimcrs 



B dimcrs 



where 



a 

IV'b) = Vs\s) + ^Va\ta) , 



(8) 

(9a) 
(9b) 



with \ta) being the three components of the triplets. This 
wave function can describe the phases that do not break 
the translational symmetry. Since we have two (i.e. A 
and B) dimers in the unit cell, the entire wave function 
l^f) is translationally invariant even when the wave func- 
tions of the two dimers are different. Certainly, this wave 
function cannot describe the plateaus except for the one 
at 1/2; for other values one would need to take a larger 
unit cell. The variational parameters u and v are then 
determined by minimizing the energy 



E = 



(*|H|*) 



(10) 



The minimization is performed numerically, except for 
some simple cases when we could find analytical solu- 
tions. 



B. Auxiliary boson formalism for the Hamiltonian 



I l'\'))/y/2, respectively. This definition is different from 
the one used in Ref. [32| by an additional phase factor 
— i which ensures that the new bosons are time- reversal 
invariant^. Furthermore, in order for the four states on 
each dimer to be faithfully represented, the number of 
bosons per dimer is constrained: 



s^s 



(11) 



The components of the spin operator at bond j are then 
given as 

= -\ dl^'^ " 4^..) - ^'o.^pjp/,, .(12b) 

The intra-dimer part of the Heisenberg Hamiltonian T-L 
[Eq. (HI)] in the bond representation reads 



H 



J 



"X ^ ^ Z ^ ^ ^a,J^a,j ' 



(13) 



J J a=x,y,z 

while the intra-dimer DM-intcraction reads 
D 



D 



■E(4 

jeB 



(14) 



We can derive similar expressions for all the terms in the 
Hamiltonian. 

In the presence of magnetic field along the z axis, it 
is more convenient to use the triplet bosons t\, ~ it],, 

and tl that create | ft), (I U) + I it))/V2, and | U) 
that are the eigenstates of the z-component of the spin 
operator. The spin operators then read: 



V2 



^/2 



V2 



^/2 



tl A 



tl 



± 



(15a) 
(15b) 
(15c) 



where the upper sign is for I = 1 spin and the lower sign 
for the I = 2 spin in the dimer, as denoted in Fig. [TJ 



In order to find the excitation spectrum, we introduce, 
in the spirit of Sachdev and Bliatl^, auxiliary bosons 
which create the singlet and the triplet on each bond; the 
operator s'f creates the singlet state (| t|) — | |t))/%/2, 
while the operators tj,, t]^, and t] create the triplet states 

*(l tt) - I U))/V2, (I tt) + I m/V2, and U) + 



C. Bond wave method 

After rewriting the Hamiltonian in terms of the bond 
operators using Eq. (|12p or (|15p , we perform a bond- wave 
approximation which is a natural extension of the usual 
spin wave theory. In the bond-wave approximation we 
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extend the number of bosons per dimer from 1 to M. so 
that the constraint (|lip now reads 



s 



E 



M 



(16) 



The variational approach mentioned in Sec. lIII Al is anal- 
ogous to finding the classical {S — ?> oo) ground state for 
spin models; M — > oo is the classical solution where 
the quantum fluctuations between the dimers are ne- 
glected. To see this, it is convenient to rotate the 'quan- 
tization axis'. Then, the variational solution ItJja) in 
Eq. ([9]) can be written in terms of the 'rotated' bosons as 
IV'a) = s^|0), where the = UgS^ +J2a '^at'l- Similarly, 
we have |V'b) = Sg|0) with Sg defined using Eq. (|9b)) . In 
the case of general M, we promote the above expres- 
sions to \iPa) = (Si)*^|0> and I^-b) = (4)^|0)' which 
are direct analogues of the Bloch coherent states for the 
spin-M/2 system. In the classical- limit Af — oo, the co- 
herent state |'0a,b) may be thought of as the condensate 
of SA,B- We also rotate the remaining bosons into ia so 
that they obey the usual commutation relations and the 
local constraint Eq. p^ . Accordingly, the expressions of 
the spin operators ([T5|) get modified. 

To consider the small 'transverse' fluctuations around 
the classical solution, we solve the constraint explicitly 
for s and treat ts as the Holstein-Primakoff bosons. Using 
the formal expansion (valid to order shown) 



= s 



(17) 



we perform a 1/Af-expansion in the spin operators and 
subsequently in the Hamiltonian. Then the Hamilton 
operator can be written as 



(18) 



where H^^^ = Ec is just the variational energy and "H*-^-* 
is a collection of the terms that are linear in 



*k ^ ( ^a;,A,k' ^y,A,k' *z,A,k' *a;,B,k' *y,B,k' *z,B,k ) (1^) 



and in the similarly defined t\^. As is expected from the 
variational nature, it turns out that "H^^-* is identically 
equal to for the variational solution. The quadratic 
part "H^^-* is of the form 



2 ^ \ t . 

ksBZ \ -k 




.(20) 



The -^(2) 

can be diagonalized by the method described 
in Appendix, Sec. [XJ and we find three excitations (one 
for each ta boson) per dimer. 

IV. BOND-WAVE SPECTRUM IN ZERO FIELD 

We start with the discussion of the zero-field excita- 
tion spectra in the low-symmetry (finite D) case. Early 
neutron scattering results^ indicated that the spectrum 
consists of essentially dispersionless (localized) single- 
triplet branch and other multi-triplet ones. This is 
the consequence of the orthogonal dimer structure, and 
triplets get dispersion in the 6'^ order of the pertur- 
bation expansion2i2£ in J'/ J. Later, higher-resolution 
neutron scattering experiments^^ revealed that the first 
triplet excitation actually splits into 3 subbands with 
well-defined dispersions. The splitting indicates the pres- 
ence of anisotropics. In the following, we will calculate 
these spectra starting from the bond-wave theory. 

The variational wave function that minimizes the en- 
ergy ([T0| in zero magnetic field takes the following form: 



I^'a) « I«) + ^ + = l-*^) + Mty) , (21a) 
w 

IVb) « I^) + ^ - = I*) - ^1*-) ' (21b) 



with 



w = £^^ = -^ + OiDyj^) . (22) 



J + V J2 + D2 2 J 
The corresponding energy is given by 



E 



Zi \Vn 



J 
2 



J2 



(23) 



This wave function is time-reversal invariant and it docs 
not break any of the symmetries of the 2?2cij the plane 
group of the Hamiltonian. We denote this phase by 

To get the excitation spectrum following the recipe 
outlined in Sec. IIII CI we rotate the states on each bond 
of type A as 



f ^1 \ 



■^A 
^.A 

\ ^Ia 



1 





w 






1 \ 

'■x.A 











1 










w 





1 
























l) 





with an analogous rotation on bonds B, so that the vari- 
ational wave functions in Eqs. pip are given as IV'a) = 
~s\\0) and IV'b) = 4|0)- Next, we condense the s]^ and Sg 
singlets. The expression of the bond-wave Hamiltonian 
is complicated for arbitrary point in the Brillouin-zone, 
except at the T point, where it assumes the following 
form: 
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1 

+ 2 



1 

+ 2 



x.B 



-2D\ 



'y,A 



x,B 
^x.A 

V Lb J 









-2D 



-2D\ 
W 



1)2 



-2D\ 









0^ 






























'■x.A 








r2 J 







-2D'^^ 
-2D'^ 



\ 






( Lb \ 




La 




La 




J 


\ Lb J 



(25) 



with 



n = 



2{D'^ + wD\^ + J'lv^) 



For simplicity, we have introduced the quantity 



'II 



II, 



D' 



lis ' 



(26) 
(27) 

(28) 



for the in-plane components of the inter-dimer DM in- 
teraction as only this combination enters the varia- 
tional ground state energy of the translationally invariant 
dimer-product wave function ([8]) and excitations. Actu- 
ally, the matrices in the Hamiltonian can be reduced 
to 2 X 2 ones by using the symmetries of the ^4 point 
group (sec Eg. IB3p . 

Following Appendix I A H we diagonalize the Hamilto- 
nian (pS)) and get the excitation energies: 



Wl,2 ~ f7 , 



= ^/n{n±2E) , 



(29a) 
(29b) 

(29c) 



We note that w J ± 2D^ and w, 



J ± 2D'^^ for 



small values of / J, thus the pairs of excitations wf 
and W4 (a = ±) arc essentially indistinguishable. Fur- 
thermore, the splitting between the two branches uj~ and 
UJ+ {n ^ 3, 4) at the F point is AD'^ + 0{D'^_^ /J ), which 
will be used to estimate the value of -D^ in Sec. IVIII Al 
A similar calculation at the q = (tt, tt) point gives 



J 



J 



J2 



^ 2 



2V J2 + D2 



(30) 
(31) 



and splittings that are quadratic in DM interactions: 

J2 _ 
2V J2 + L'2 4J2 ■ 



UJ2 — Wi 



J 
2" 



(32) 



Let us note that for larger values of | D'^ | the dispersion 
becomes comparable to the gap, and new phases appear. 
The branches become gapless when \/ ,P + D"^ — 
T^{D'^ + wD\^ + J'w2), while = for 4L»^ = 

=fV J2 4- £)2^ Assuming that Dy is absent and keeping 
only the leading term in D/J, we get that the phase 
ZiYD2d\ is stable for 



J 

4 



8J ^4 



El 

8J2 



(2J' - J) 



(33) 



in the zero field. Beyond these boundaries Z2 twofold 
degenerate phases with the symmetry group C2V (when 
i^J^ for D'^ > 0) or ^4 (when ^ for < 0) is 
realized (sec also Fig. (5]). Wc will discuss these Z2 phases 
in more detail in a separate publication. 




(0,0) 



M (Ji,0) 



(0,0) 



FIG. 3. (color online) Dispersions of the quasi-triplet excita- 
tions in zero magnetic field: cu^^ (bottom), a;i,2 (middle), and 
oj^^ (top). We have chosen D — D'^ = O.IJ and J' — 0.6 J. 



After a first order expansion with respect to D/J^ 
D'^l J , Z?!! ^/J, and _D|| it is possible to solve the 

eigenvalue problem analytically. In this case wc get three 
two-fold degenerate branches: a dispersionless one with 
eigenvalue J and two branches with 



j2 ± 



(34) 
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where 



-^'^ 2D' 



(1 - cos go cosgfa) 



+16DY cos' 



qa 2 96 
— COS — 

2 2 



1/2 



(35) 



The dispersion of the quasi triplet excitations is shown 
in Fig. [1 

The dispersion in zero field has also been considered 
by Cheng et al. in Ref. [1^1, where they used a different 
approach; by suitable rotation of the spin operators they 
removed the D and arrived at an effective Hamiltonian, 
where they carried out a first order perturbation expan- 
sion to get the dispersions of the effective triplets. Even 
though they considered a unit cell that has C4 symmetry, 
our dispersion agrees with their result, up to ambiguity 
in the sign in front of the D'^^ ^ in Eq. ((35|). Furthermore, 
they extended their analysis by exact diagonalization cal- 
culations of the spectra. 
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In this section, we are going to consider the variational 
ground-state phase diagram in the presence of an exter- 
nal field along the z-axis. The full Hamiltonian is now 
invariant under the magnetic group S^ + Qaxz^S^ which 
is isomorphic to 'D2d- For clarity of the argument, we 
investigate the high-symmetry case (where D = and 
gs = 0) and the low-symmetry case (which realizes in the 
low-temperature phase of SrCu2 (603)2) separately. 



FIG. 4. (color online) (a) Phase diagram in the h^-D'j^ plane 
for D — and J'/ J = 0.6. DS is the dimer singlet phase that 
remains a variational ground state even for finite values of D(l • 
The spin configurations in the different phases are shown in 
Fig.O 'm=l/2 plat.' denotes the half-magnetization plateau 
phase, with a singlet and a magnetized triplet in each unit cell, 
(b) The magnetization curves for different values of \D'j_\ as 
a function of the magnetic field. 



A. High symmetry case 

When the space group is I A/ mem (which is relevant in 
the high temperature phase T > Tg) and the symmetry 
group of the two-dimer unit cell is V^h, the intra-dimer 
DM interaction D is absent and only the inter-dimer D'j^ 
DM-interaction is finite. With this type of anisotropy 
the z component of the spin is a conserved quantity and 
this greatly simplifies the form of the variational ground 
states and of the bond-wave Hamiltonian. 

Numerically minimizing the variational energy 
in the presence of a magnetic field 
along the z direction, we have found three gapped 
phases (see Fig. |4|): the dimer-singlet (DS), the one-half 
magnetization plateau, and the fully polarized phase. 
Furthermore, there are four gapless phases associated 
with the symmetry breaking of the continuous 0{2) 
symmetry: the Neel, the 0(2)[C4], the 0(2)[54], and 
the 0(2) X Z2 phase. In these 0{2) phases, the 0(2) 
rotational symmetry in the xy plane perpendicular to 
the field is spontaneously broken and they are the 
consequence of the S^' being a good quantum number. 
Below we will consider the different phases and their 
excitations in more detail. 



1. Dimer-singlet phase 

As we mentioned earlier, the exact ground state of 
the SU(2) symmetric Shastry-Sutherland model is the 
product of singlets on dimers: \iPa) = IV'b) = \s) for 
< J' < 0.68 J, as shown in Ref. [3]. In the variational 
approach this ground states turns out to be stable for 
finite values of D'^ and magnetic fields h < h^, where the 
critical field is given by 

h, = ^J^~4\D'^\J . (36) 

The ground state energy is coming purely from the ex- 
change within a dimer: 

Ebs=-Y' ^^^^ 
all the other bond energies are identically 0. 

2. The 0(2) phases 

Between the dimer singlet and the one-half magnetiza- 
tion plateau we find 0{2) symmetric phases (see Fig. S]) 
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0(2)[C4] 



0(2)xZ2 




0(2)[S4] 



Zi[D2d] 



FIG. 5. (color online) Schematic representation of the 0(2) 
symmetric phases in the high symmetry [D = 0) case with 
magnetic field perpendicular to the CuBOs layer (i.e. h\\z). 
We have plotted the expectation values of the spin component 
in the xy plane. Since in this case is a conserved, the spins 
can be arbitrarily rotated by a global 0(2) rotation in the 
plane. The blue (open arrow) and red (solid arrow) represent 
inequivalent spins. 



when we apply the field perpendicular to the plane. In 
these phases the magnetization increases continuously 
between and 1/2 per dimcr (or m~ ~ 1 per unit 
cell). Since the is a conserved quantity, the Hamilto- 
nian does not break the 0{2) symmetry of the rotations 
around the z-axis. This symmetry is spontaneously bro- 
ken in the 0(2) phase. From numerical minimization we 
found that the wave function can be written as 

IVa) oc \s) + ue^^lti) + de^"^|ii) , (38a) 
IV-b) oc \s) ± iue''P\ti) T ide-''P\ti) , (38b) 

with the upper sign for D'^ > and the lower sign for 
D'^ < 0. This wave function is continuously connected 
to the dimer-singlet phase, as setting u = and d = 
we get back to the product of singlets. In this phase 
the expectation values for all the spins are equal 
{{S^} oc |wp — |(ip), and the spin components in the 
xy plane along the intcr-dimcr bonds arc perpendicu- 
lar to each other in such a way that going around on a 
void square the spins also make a full turn (as shown in 
Fig. [5]) - in other words, the spin configurations are in- 
variant with respect to either the C4- or the S'4-rotation. 
It is the sense of the rotation that makes the difference 
between the two different cases of the sign of Z?^. To 
make a clear distinction, wc use the symmetry groups 



which leave the variational ground state invariant to la- 
bel these two phases as 0(2)[C4] and 0{2)[S4]. We found 
that they are realized for positive and negative values of 
D'j_, respectively. 

The expectation value of the Hamiltonian with the 
wave functions ^38a\i and (|38bp reads 



Eo(2) — 



4{u - df 



d^-3J , 2{u^ - dy 

(w2 + d^ + If 

J -Pi I- ; ,(39) 



and the minimization gives a set of polynomial equations 
that needs to be solved numerically. Close to the phase 
boundary to the dimer singlet phase given by Eq. (|36p . 
we can expand in Sh ~ — he- In the lowest order in 6h 



2^AJJ'hl + J^~h^^ 
(J - hc)y/2h^ 



Sh 



(40a) 
(40b) 



2^AJJ'hl + J^-ht 
The magnetization below /ic is 0, and above he grows as 

2J\D'^^\-AD'^^ + J'{J-A\D'^\) ^ ' 

(41) 

In the absence of the magnetic field the modulus of 
the amplitudes u and d of the two triplet components 
become equal, and writing v/y/2 = u = ~d the wave 
function simplics to 



\ti)) , (42a) 



|^B> « \s) ± —I {e^\t^) + e-'^\t-^) , (42b) 

where, as we noted, the sign is determined by that of 
D'j^ = ±\D'^\. The minimum of the energy Eq. in 
this case is achieved for 



I4\D'^\-J 



with 



E 



J 

2 



2\D\ 



8\D\ 



(43) 



(44) 



From the analysis above it turns out that the 0(2) phases 
arc realized for \D'j_\ > J/4 in zero field [this is consistent 
with Eq. (p6t]. 



3. The 0(2) x Z2 phase 

The 0(2) phasc(s) and the one-half magnetization 
plateau phase are connected via two continuous phase 
transitions. The intermediate phase exhibits both the 
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Z2-synimetry breaking of the plateau phase (the inequiv- 
alence of the z-component of the magnetization on the A- 
and the B dimers), and the 0{2) symmetry breaking of 
the 0{2) phase, as shown in Fig. [SJ As we approach the 
boundary of the one-half magnetization plateau the com- 
ponent of the spins perpendicular to the field decreases, 
and eventually vanishes at the phase boundary. Though 
we do not break the translational symmetry, the fact that 
the magnetization along the field is not equal on the A- 
and the B dimer (a discrete symmetry is broken), and 
that at the same time we break a continuous symmetry 
of the 0{2) type, we may call this phase a supersolidi^i^i 
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1 ' 1 ' 1 ' 1 
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,1,1,1 


0.2 






0.1 



-0.1 


Zi[D2d] 




-0.2 
-0.3 


1,1,1,1,1 


1,1,1, 



0.5 



0.4 
0.3 



0.2 



0.1 



:(b) 


















1,1,1,1,1,1,1, 



0.2 0.4 0.6 0.8 



1.2 1.4 1.6 



FIG. 6. (color online) (a) Phase diagram in h^-D'^ plane the 
for J' / J = 0.6 and D/J = 0.1 (the low-symmetry case). In 
comparison to the D = Q case in Fig. U in a large region 
of the phase space the dimer-singlet and the 0(2)[S'4] essen- 
tially merged to create the Zi [D2d\ phase, a small part of the 
0(2)[S4] phase become a twofold degenerate ^2(54], and the 
0(2)[C4] merged with the m = 1/2 magnetization plateau 
phase into the Z2[C2v\ phase, (b) Magnetization curves for a 
few selected values of D\ . 



B. Low— symmetry case 

At low temperatures (specifically, T < Ts — 395 [K] 
in SrCu2(B03)2), the symmetry of two-dimer unit cell is 
I?2d- As noted earlier, the lowering of the symmetry al- 
lows for a finite value of the intra-dimer DM-interactions 
D. From numerical minimization, we mapped out the 
phase diagram, and we found a ground state that does 
not break any of the symmetries of the Hamiltonian in 



the low field region of the experimentally relevant param- 
eters (we call the phase Zi['D2d]) (see Fig. Addition- 
ally, we found many two-fold degenerate Z2 phases that 
we will describe in more details in a separate paper. 

The total is not a good quantum number any more 
and the continuous symmetry of the 0(2) phases gets 
reduced to discrete symmetries, and in the absence of 
continuous symmetry all phases become gapped. This 
symmetry reduction is seen in the expectation values of 
the energy; the inclusion of the intra-dimer DM D and 
the staggered g-tensor gs introduces anisotropy terms to 
Eq. (133 



E 



2i[X'2d] 



0{2) 



V2D- 



u 



dp + 1 



cos If 



U2 + rf2 + 1 



COS (p 



(45) 



which determine the preferred direction of the xy com- 
ponents. Assuming D > and positive u and v, the 
DM energy on the dimers is minimal when ip = t:. The 
variational wave function in the Zi \D2d] is given by 

I^a) « \s) - u\h) - d\ti) , (46a) 
IV'b) oc \s) ^ iu\ti) - id\ti) , (46b) 

and the expectation values of spin components are shown 
in Fig. [5l This phase is adiabatically connected to the 
dimer-singlet product phase (when it — >■ and d 0) 
and in general the intra-dimer DM interaction D mixes 
the triplet components to the singlet, as expected from 
Eq. (|T4|) . It is also a special case of the wave function of 
the 0(2)[54] phase [ Eqs. ([55a|) and (|38b| with lower sign] 
with the phase locked to = 0. In other words, when 
D' < Q (the experimentally relevant case) turning on 
an infinitesimal value of D removes the phase boundary 
between the dimer-singlet (DS) phase and the 0(2)[iS4] 
phase. For D' > 0, the 0(2) [C4] phase becomes frustrated 
with respect to D, and will give rise to a Z2-symmetry 
breaking. 

The conditions dEz^[D2a]/du = and dEz^[V2d]/ ^'^ = 
lead to a set of polynomial equations of high degree that 
one can solve only numerically. However, we can search 
for the solution of u and v as an expansion in D/J. For 
small values of D, both u and d are linearly proportional 
to D when hz < he [he is defined in Eq. ([36])], and we 
can expand the energy as 



E 



Zl [V2d] 



3J 

T 



+ 2J{d^ + u"^) + 4L>^(ii - d)2 



-2hz{u^ - d^) - V2D{i 
-2\/2gshz{u — v) , 



d) 



(47) 



and in the lowest order in D and g,, the minimum is 
achieved with 



1 {D-2gshz){J + hz)+^DD'^ 
271 J2 + UD'^ - hi 

1 {D + 2~gsK){J -hz) + ADD'^ 
2^/2 + UD', - hi 



(48a) 
(48b) 
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For hz = we recover Eq. ([22]) . To be more precise, the 
expansion is actually in D/{,P + 4:JD'^ - hi) = D/{hl - 
hi), and the denominator becomes at the I? = 0, = 
boundary between the dimcr-singlet and the 0(2) [^4] 
phase. The expression for the magnetization of a dimer 
is then 



VI. BOND-WAVE SPECTRUM FOR THE 
FIELD PARALLEL TO z AXIS 

In this section, we calculate the bond-wave excitation 
spectrum in the presence of an external magnetic field 
h \\ z and examine qualitatively the effect of intra and 
inter DM interactions. 



A. High symmetry case 



J {hi-hiy 



(49) 



which grows quadratically with the anisotropy. We note 
that J' enters only in the next order in the expansion. 



First, we consider the case where D = Q. In the dimer- 
singlet product state we condense the singlet states on 
both the A and B bonds, and there is no need to rotate 
the bosons [i ^t). The energy is then (up to a constant 
shift) 



n = Ej,s+Y. ^2(q) , 

qGBZ 



(50) 



where 



-H2(q) = J 4.A(q)^o,A(q) + 4,B(q)*o,B(q) 
^ 4,B(q) ( J 



•i,A(q) 
q) 



h.Ai 



^21?;^ 71 J-hz i2D'^ji 



i2D'j^ji 

Vo 



J + h 



2iD' 



'l7i 



2iD'j_"/i J + h 



( 



1,B 



tt 



(q) \ 

(q) 



1,A 



(51) 



and 



71 



qa qb 

cos — cos — 
2 2 



(52) 



The Hamiltonian matrix is of the form of Eq. ()A6p , and 
can be diagonalized following the procedure outlined in 
the Appendix IA 21 The operators in the momentum space 
are defined by the t^(k) = N^^^'^ e*'^ '"jt], where rj is 
the position of the j-th spin. 

Actually, we can introduce the following combinations 



il,±(q) = ^ [4,A(q)T*its(q) 



(53a) 
(53b) 



together with the corresponding annihilation operators, 
so that the original 4 by 4 matrix in Eq. (|?T|) decomposes 
into two 2 by 2 problems, with the Hamiltonians 

Hi''(q) = {J-K.± i?l)iX±(q)ti,±(q) 

±2^171 [il±(q)tL(q) + 4.±(q)%,±(q) 

+ (J + /i, ±i?l)t|^(q)fj_^(q) . (54) 



The bond-wave spectrum consists of six modes: 
twofold degenerate nondispersing excitations with 
w(q) = J (denoted by Tq'° in Fig.[7|) and four dispersing 
modes: 



-t = aI ,P + AJD'j^ cos — cos — ± h^ (55) 



that come from "H^ ' (q) (blue solid lines denoted by T^i 
in Fig. [7]) and 



4:JU , COS — COS — ± h, 
^22 



J2 



(56) 



from "HIl (q) that we will call the T^^ modes (shown by 
red solid lines in Fig.[7|). The dispersions have a finite gap 
in the dimer-singlet product state. Let us also mention 
that for /i^ = and Z) = we recover Eq. ((35|) . 

From the equations above, the gap closes at q = 
when the magnetic field reaches he defined by Eq. 
and we enter into the 0{2) phases. Unfortunately, the ex- 
plicit form of the variational wave function and the bond- 
wave Hamiltonian in the 0(2) phases is too complicated, 
thus here we discuss the numerical solution only. We just 
mention that the combinations t\ _|_(q) and tj _|_(q) intro- 
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duced in Eqs. (|53)) are decoupling the bond-wave Hamil- 
tonian in the 0{2) case as well. When D'j_ > 0, the 
closing of the gap leads to a Goldstone mode — the con- 
sequence of the continuous symmetry breaking — that 
appears as a continuation of the mode, and the 

condensation of a linear combination of the t[ _{q = 0) 
and tj _(q = 0) bosons results in the 0(2)[C4] phase de- 
scribed by the (|38p wave functions with the upper sign 
(see also Fig. U]). For Z?^ < 0, on the other hand, the 
Goldstone mode evolves from the cj+._ mode, leading to 
the 0(2)[54] phase (Fig.H])). As we can see in Fig. [3 the 
lowest gapped mode for q = in the dimer-singlet phase 
remains gapless while the 0(2) symmetry is broken and 
this is the case until the half magnetization plateau. 

Also, from the q-dependent excitation spectrum we 
learn that the spectrum may become gapless not only 
at the q = 0, but also at some other q values in the Bril- 
louin zone, thus announcing a helical instability of the 
0{2) phases. In Fig. [7] we indicate the boundary of this 
instability (the hatched region) that we have obtained 
from the numerical calculations of the spectra. 

The strength of the magnetic probe response is de- 
termined by the structure factor S'""(q, w). In particu- 
lar, the S^^{q = 0,aj) and S'^^(q=0,a;) determines the 
strength of the ESR lines in first approximation, when 
the static magnetic field is along the z axes. The struc- 
ture factor is given by 



5""(q,c^) « J2 |(/l^q|0)r Si"" -Ef+ ^o) , (57) 



where |0) is the ground state (in our case the |^') varia- 
tional wave function), by / we denote the excited states, 
and Eq and Ef are the energies of the respective states. 

As a first step, it is instructive to look at 
the w-integrated (static) structure factor, ^""(q) = 
/ ^^^""(q, w), which is actually the sum of the (positive) 
matrix elements, and that is equal to (>I'|S'"q5q|^'). In 
the (pure) dimer-singlet ground state, limq_>.o "SqlO) — > 0, 
so we expect to see no response in ESR experiments, un- 
less there are anisotropics which mix the triplet compo- 
nents with the singlet. 

In the 0{2) phase (discussed in Sec. IV Al 2). the static 
structure factor is 



5""(q = 0) = 



4JJ'/l2 + J4 _ M 



Sh 



(58) 
(59) 



for a = a;, y, and z. Examining the individual matrix 
elements, it turns out that the matrix elements for the 
S"^^ and 5*^^ are all vanishing except for the Tg line. 
On the other hand, the matrix elements for the are 
nonzero for the Tf and Tf lines. Since the ESR line 
width is proportional with S^^ and S*^^ when the field is 
along the ^-direction, we expect a strong signal for the 
T° fine. 
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FIG. 7. (color online) Bond-wave spectrum at the F point 
when we keep the singlets and the S'' = \ triplets only. The 
magnetic field is parallel to z-axis, J' = 0.6J, and D'^ = O.IJ. 
The transition into the 0(2) phase happens at hz ~ 0.77J, 
and into the plateau phase at « l.lj. For 0.77 < h./J < 
1.03 one of the excitations becomes identically zero — this is 
the Goldstone mode of the 0(2)[iS4] phase. The Goldstone 
mode of the 0(2) x Z2 phase is zero for 1.03 < h^/J < 1.1. 
The 2-dimer variational solution is unstable in the gray region 
- the dispersion goes to at some wave vector away from the 
F-point at hz = 0.96 and 1.08. The filled area above the 
dispersion line shows the strength of the spin structure factor 
S^^+S^y. The dotted line is the approximation from Ref. [H]. 



B. Low symmetry case 

Firstly let us discuss the case D'j^ = shown in Fig. 
[HI At low fields, the spectrum looks like the usual one- 
triplet excitation Zeeman-split by the magnetic field: 
we see three branches that arc two-fold degenerate as 
we have two dimcrs in a unit cell. Without any kind 
of anisotropics, these excitations would correspond to 
the pure one-triplet excitations. However in the low- 
symmetry case the intra-dimer DM coupling D mixes the 
singlets with these excitations. From the zero-field equa- 
tions Eqs. we get that the splitting is of the order of 
D- for small values of D/J, which is much smaller than 
the linear splitting caused by D'j^. On the other hand, 
the effect of D is much more pronounced at higher fields, 
where the gap becomes small and the singlet-triplet mix- 
ing is enhanced. Instead of the Goldstone mode, the 
anisotropy induces a "level repulsion", and we get a fi- 
nite gap that is roughly proportional to ^/dJJ, that is 
consistent with the usual form of the anisotropy gap. We 
note that the "level repulsion" happens only to one of the 
two almost degenerate branches that come down with the 
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applied field, and it depenctS, .officially on the symmetry 
of those state as has been noted by Miyahara and Mila in 
Rcf. (fssh. where they considered the dispersion of a sin- 
gle triplet bond moving in the singlet background by the 
standard perturbation theory. As we increase the field, 
the gap closes for the Tf level at the phase boundary to 
the Z2 phase with C2V symmetry (see Fig.[S]). 
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FIG. 8. (color online) Excitation spectrum in magnetic field 
parallel to z-axis when _D(l = 0, J'/J = 0.6. The instabilities 
toward helical states are at /iz/J = 0.9421 and 1.1002, while 
the k = instability into the Zi phase is at 0.9515. 



For finite inter- and intra-dimer DM interactions, we 
observe both the zero-field splitting and the anti level- 
crossing around the critical field. In Section IIVI we 
presented a detailed calculation for the zero-field disper- 
sion and estimated this splitting in the first order of DM 
interactions to be 41?^. This is in excellent agreement 
with the findings of Cheng et al [s^. The two low ly- 
ing modes in Fig. [9] and ITO)) curve differently in the 
0(2) phase, only one of them crosses the ground state 
and while the other is gapped. This can be explained by 
that only one (the Tf ) low lying excitation is coupled to 
D, and the gap is proportional to it. [11]. In the case of 
the bond-wave calculation, the gap opens as \/D as the 
effect of quantum fluctuations (see Appendix |B] for more 
details). Flipping the sign of the inter-dimer DM cou- 
pling D'j^ changes the lowest-lying mode (compare Fig. 
|9]and [TO)) . The singlet-triplet mixing is different accord- 
ing to the symmetry of the lowest-lying mode and the 
anti level-crossing occurs only for D'j^ < 0. 

The Eq. ([55)1 for the static structure factor is valid also 
for finite D values with the variational parameters u and 
V in the wave function (|46p now being obtained by the 
minimization of Eq. (j45|) . In the limit of small D/ J, we 
can use the values that are given by Eqs. pS)) . and at 
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FIG. 9. (color online) Excitation spectrum in magnetic field 
parallel to z-axis in the low-symetry case, for D'^ > 0, D/ J = 
0.1, and J'/J = 0.6. The notations are the same as in Fig. [T] 




FIG. 10. (color online) Excitation spectrum in magnetic field 
parallel to z-axis for D'j^ < 0. (J' = 0.6J) 



q = we get in lowest order in D/J 

{J + AD'^_f + hl 



(60) 



4 [j-^ +UD'^-hiy 

for small values of field (the apparent singularity at the 
critical field ^ = J'^ — 4|Z?^| J is an artifact of 
the expansion). Similarly to the D = case, the weight 
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of the spin correlation function S'^^(q = 0,uj) in the 
^1 [2-'2d]-phase is concentrated on the single Tq line. As 
we enter the Z2[C2v] phase, the ^^^(q ~ 0,aj) is split 
between the Tq and Tq lines. The strength of the indi- 
vidual lines is shown in Figs. HI |8l and [10] by the filled 
region above the Tq and Tq lines. 



VII. PHASE DIAGRAM AND EXCITATION 
SPECTRUM FOR IN-PLANE MAGNETIC FIELD 

In this section we consider the case when the magnetic 
field is applied in the CuBOa plane, parallel to the bond 
A and perpendicular to bond B. This direction is denoted 
by X in Fig. [T](we note that this direction is different from 
the crystallographic a-axis). This choice makes the two 
dimers inequivalent and breaks the rotational symmetry 
S4. This direction of field lowers the 'D2V symmetry of the 
unit cell in the low-temperature phase to the magnetic 
group {E,<7yz} -\- 9C2(z) X {E,ayz} that is isomorphic 
to C2V In the following, we show the phase diagrams 
for finite values of D and give a short discussion of the 
phases that appear. We also show the ESR spectrum in 
the end of this section. 



A. Phase diagram 



(Si X S2) = —Vx/il + vl). This is the so called p-type 
ncmatic state<^i^ The parameter Vx does not depend 
on the magnetic filed, and minimizing the energy (p^ 
we find that = D/2J. 
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FIG. 11. (color online) Phase diagram for the magnetic field 
in the x-direction. J' /J — 0.6 and D/J — 0.1. The spin 
configurations in the difi'erent phases are shown in Fig. [12] 



The numerically obtained phase diagram as a function 
of hx and D'j^ for a selected value of D is shown in Fig. 
[TT] We denote the ground state by Zi[C2v], which has the 
full C2V symmetry of the Hamiltonian and the variational 
wave function is of the following form: 

IiPa) (X \s) - Vy\ty) - iuz\tz) , (61a) 

\4>b) \s) + Vx\tx) , (61b) 

with the energy expectation value: 

-E'Zi[C2„] 

J(l - vl) + 2Dvx J + 2hxu,vy + Du, 

2(1 + 1.2) 1 + ul+vl ■ ^ > 

The intra-dimer DM interaction D prefers states with 
dipole expectation values that are perpendicular to the 
vector D. As a consequence, the magnetic field induces 
the moment only on dimers where the direction of D is 
perpendicular to the field. In our case h\\x^ and only 
the dimer A develops a finite magnetization: the spin 
components along the magnetic field, S% i = Sj^ 2 '^^11 
as the components perpendicular to the plane, ^ = 
2 become finite, as shown in Fig. [T^] Increasing 
the magnetic field, the expectation value of increases 
smoothly up to — 1/2 [sec Fig. [T3] (a)]. 

The wave-function {iPb), on the other hand, is time- 
reversal invariant, where the expectation value of any 
spin component is zero: (81^2) = 0. However, it breaks 
the rotational symmetry, as the vector chirality is finite: 




FIG. 12. (color online) Schematic representation of the spin 
configurations in the phases Zi [C2v], Z2 [Ci^'] and Z2 [Cs'''] that 
appear when the field is along the x direction. The darker and 
lighter arrow represent the two degenerate states. 



On the phase diagram, we found two additional 
phases, and they are both two-fold degenerate Z2-phases. 
For sufficiently large negative values of D'j^, the phase 
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Zi\p,s \ is realized with the wave function 

I'/'a) « |s) - {vy ± iuj^) \ty) - [iu^ T Wz) \iz) , (63a) 

IV'b) OC |.s) + (Vx ± iUa;) \tx) ■ (63b) 

The magnetization on the A dimer consist of a uniform 
part S"^ 1 = 5"! 2 ^^'^ ^ staggered part -S*^ ^ = "^a-, ^^'^ 
S^-^ = —8^2 shown in Fig. [T2). While in phase 
Zi[C2v] there were no dipole components on dimer B, 
here the expectation value of staggered magnetization 
to|j St is non-zero, ^ ~ ^^^b 2- '^^'^ magnetization 
pattern is invariant under ayz- 

At large enough positive D'j^, we reach the phase 
■^2[Cs'^''] where the ground state has the following form 

= \s) ± iux \tx) - Vy \ty) - iuz \tz) , (64a) 

|*>B = + \tx) ± lUy \ty) ± V, \t,) . (64b) 

In this case, the spin expectation values are, as is shown 
in Fig. 1121 invariant under the reflexion Quxz ■ 

We note that in the limit of Z? the Zi[C2„] -phase 
is continuously connected to the dimer-singlet phase that 
is realized for 

hx < \J.P - IQD'^^ (65) 

when D = Q. It is also continuously connected to one 
of the ground states of the twofold-degenerate to = 1/2 
plateau phase (where the singlets are located on the B 
bonds). 

B. ESR spectra 

In the following we discuss the effect of the DM compo- 
nents on the ESR spectrum for the magnetic filed parallel 
to a;-axis. We remind the reader that in the absence of 
the anisotropy (i.e. DM interactions) the dimer singlet is 
the ground state for low fields hx < J and that the exci- 
tations arc the pure, Zeeman-split triplets with energies 
J — hx, J, and J + hx each of which is twofold degenerate 
corresponding to the two dimers in the unit cell. 

In Fig. [131 show the calculated ESR spectrum for 
D = O.IJ and D'j^ = — O.IJ. In the absence of the field 
we observe the zero field splitting 4Z?^ that we discussed 
in Sec. IIVI Now, unlike the case of the field along the z- 
direction, the spectrum consists of three pairs of almost 
degenerate levels (note that in the absence of D each 
pair is exactly degenerate in the dimer singlet phase), 
and only at higher fields near the the phase transition 
the lines split. When D'j_ is large enough, the gap closes 

at the boundary to the Z2[ci^^]-phase. For = O.IJ 
the spectrum looks essentially the same. 

For D'j^ = (i.e. only D is present), the zero- field 
splitting disappears and as we approach the critical field 
the two-fold degenerate triplet branch splits. The difi^er- 
ent behaviors of the two low-lying excitations around the 
critical field is due to their different singlet-triplet mix- 
ing. As has been discussed previously when D is finite 



and D'j^ is zero we are in phase Zi[C2v] (see the phase di- 
agram Fig. [TT|) . Increasing the magnetic field from zero, 
the value of the parameter Uz in the ground state wave 
function \iPa) (Eq. (HU) increases continuously devel- 
oping a finite magnetization at dimer A, while the 
magnetization of dimer B remains zero. 




FIG. 13. (color online) Excitation spectrum in magnetic field 
parallel to x-axis. (J' = 0.6J) 



VIII. COMPARISON TO EXPERIMENTAL 
SPECTRA 

The ESR spectra have been considered previously 
by perturbation^^ (for D = 0) and by exact 
diagonalizations^. In our approach, it is straightfor- 
ward to take all the anisotropics which are relevant in 
experiments into account and below we consider the ESR 
spectra in a more realistic setting to test our theoretical 
framework. 

Various attempt have been made to determine the val- 
ues of the different terms in the Hamiltonian. For com- 
pleteness, we write down the Hamiltonian in its full form: 

n = j^s.Sj + J' J2 ^^^j 

n.n. n.n.n. 

+ ^D,j (S, X S,) -f D',, • (S, X S,) 

n.n n.n.n 

— ^ Qzhz {S\i + + •S'^i + ^m) 

u.c. 

+ J2 9shz {SI, ~ + SI, - SI2) . (66) 

u.c. 

Using the estimations of Cepas et al. (Ref.[l^) for 
the value of the inter-dimer DM component (D'^j J = 
—0.02), Kodama et al determined the intra-dimer DM 
component via fitting exact diagonalization data and 
obtained^^ D/J = 0.034. Fitting the neutron scattering 
data at zero field^, Cheng et al. predicted^ Z?^ = 0.18 
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meV and D[| ^ + J'D/2J = 0.07 mcV. Using ah ini- 
tio LSDA+U calculation, Mazurenko et al. estimated^ 
D = 0.35 mcV for intra- and D'^ = 0.1 mcV, D'^^ ^ = 0.06 

mcV and Z^ii — 0.04 mcV for intcr-dimcr components. 

The values of the (^-tensor anisotropics were estimated 
from ESR and NMR measurements: gx = gy = 2.05 and 
g^ = 2.28 in Ref. [l^], and g^ = 0.023 from the tilt angle 
of the electric field gradient in Ref. . 

Information on the excitation spectra are avail- 
able from ESR^i^ Raman^ii^ far-infrared (FIR) 
spectroscopy r^i^ and neutron scatteringi2i22ii2i measure- 
ments. We will mainly compare our spectra with the 
ESR, far-infrared, and neutron-scattering measurements. 
Our lines T^^ correspond to Toj,(±) in the FIR spectra 
in Ref. [H and to d in the ESR spectra in Ref. [H], 
the lines T^^ correspond to Tom(±) in [i^ and to O2 
Ref. [H, and the lines T° and T° are Top,™(0) in HI]. 



A. Quantitative comparison to experiments at zero 

field 

In particular, high resolution ESR measumerents 
of Nojiri et al^ sees the two triplet excitations 
at 679 ± 2GHz and 764 ± 2GHz. The FIR mea- 
surements of Room et aZ4^ observed three triplet 
modes at 22.72±0.05cm-i («681GHz), 24.11±0.05cm-i 
(w723GHz), and at 25.51±0.05cm-i (w765GHz). The 
origin of the signal at 24.11cm~^ is the A 5*2 = triplet 
excitation that is not seen in the zero field ESR spectra. 
From Eq. ([55]) we can deduce that the splitting between 
the AS*^ = 1 and AS"" = -1 triplet hues gives 4i?^ « 85 
GHz, that is D'j^ « 21 GHz. 

Furthermore, high resolution inelastic neutron scatter- 
ing measurements in zero field performed by Gaulin et 
al^ have revealed that the dispersion above the gap 
consist of three distinct branches of triplet excitations. 
The splitting observed there has been fitted by Cheng 
ali^ to yield the result which is identical to our Eq. ([55]) . 
From these dispersions, the splitting between the triplets 
at q = is 17(o,o) = 4:D'j_ « 0.4meV (that is «95 GHz, 
close to the 85GHz given above), while at q = (tt, 0) it is 

r!(,,o) = V2 (2i?[| , - ^) = 0.2meV. 

We need to mention here that in our approach the 
dispersion of the triplets coming from the inter-dimer 
coupling J' is altogether missing; this is why we have 
used, in estimating D'j^, the 'bare' value J of the single- 
triplet gap, instead of using the renormalized value which 
is actually observed in ESR measurements (the above 
value may be renormalized if we go beyond the linear 
bond- wave approximation). Numerical diagonalization 
of Cheng ali^ has shown that the effect of J' is only to 
modify the dispersions in such a way that the splittings 
remain independent of the J' . 



B. Quantitative comparison of the spectra at finite 
magnetic field 

In this subsection, we try to fit the ESR spectrum of 
Nojiri et al. in Ref. j23| by using the bond-wave method 
on top of variational calculation. To obtain a quantita- 
tively good fit, we need to include the DM interactions 
as well as the y-tensor anisotropics . 




5 10 15 20 25 30 35 40 
h, (T) 

FIG. 14. (color online) Qualitative comparison of excitation 
spectrum witli the h\\c ESR spectrum shown in Fig. 4(a) in 
Nojiri et a/. [2^. In the inset (b) the dimanonds are the far- 
infrared data from Ref. [4^ , and the squares are the ESR data 
hom Ref. [13] 



In the fitting, we use the values of the anisotropy con- 
stants gz_= 2.28 estimated in Ref. [l^l, and ~ 0.023 
in Ref. 125|. The value of the intra-dimer DM coupling 
D = 0.034 J « 60 GHz is obtained in Ref. assum- 
ing that J = 85K (Ref. and [ll| , and a similar value 
{D = 1.8 cm-i = 54 GHz) is reported in Ref. jH. The 
inter-dimer coupling constant is given by Z?^ = 21 GHz, 
as determined above. For the reason described above, we 
choose J to be equal to 722 GHz, the value of the ex- 
perimentally observed gap2ii. Furthermore, we find that 
the spectrum is essentially independent of the value of J' 
inasmuch as we are in the Zi\D2d] phase, so we have cho- 
sen J' / J = 0.6 for internal consistency of the calculation. 
The calcualtcd bond-wave spectrum with the parameters 
mentioned above is shown in Fig. 1141 

We find a surprisingly good quantitative agreement 
with the high- field ESR of Nojhi et al (Ref. Jp|) and 
the FIR measurements of Room et al. (Ref. jl^). Our 
spectra reproduced not only the value of the high-field 
gap in the Tf-excatiation above 20T, but also the be- 
havior of the Tq level which follows nicely the main (i.e. 
largest- intensity) peak in the ESR spectrum [Fig. [TU^b)], 
thus clearly identifying those lines as originating from 
triplet excitations. 
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IX. SUMMARY AND CONCLUSIONS 

We have studied the description of the magnetic prop- 
erties of SrCu2(B03)2 using the Shastry-Sutherland 
model extended with anisotropics. The possible form of 
the anisotropics, like the Dzyaloshinskii-Moriya interac- 
tions [Eqs. (|4])-(l5])] and the g-tensor anisotropy [Eq. ©J 
follows from the structure and the symmetry properties 
(space group) of the material. 

We used a bond-factorized form of the variational wave 
function to study the effect of the anisotropics on the 
ground state properties in phases which are compati- 
ble with the crystallographic unit cell comprising of two 
orthogonal dimcrs in the presence of external magnetic 
field. This includes the experimentally relevant dimer- 
singlet phase in low fields below the 1/8 magnetization 
plateau and the 1/2 plateau. 

We have found that the dimer-singlet phase remains a 
good variational ground state in the so-called high tem- 
perature phase, where the only anisotropy, that takes on 
a finite but small value, is the inter-dimer interaction 
D'j^ perpendicular to the CuBOa planes, and the mag- 
netic field is perpendicular to the plane. This phase is 
surrounded by gapless phases with an 0(2) symmetry 
where, due to D^, the triplons can propagate with a 
well-defined dispersion. 

In the less symmetrical, low-temperature structure of 
SrCu2(B03)2, the finite intra-dimer DM interaction D 
appears and it gives rise to admixture of triplet compo- 
nents and the singlet states in the variational wave func- 
tion of the 'dimer-singlet' phase, without any symmetry 
breaking (we denoted this phase as Zi\D2v\)- The finite 
D not only removes the 0(2) degeneracy in the gapless 
phases, but also introduces frustration depending on the 
sign of D'^-. in the unfrustrated case there is a crossover 
from the 0(2) phase to Zi\V2d\ where there is a preferred 
direction in the 0(2)-plane set by Z?, while in the frus- 
trated case an Ising like Z2-symmctry breaking appears. 
This is a similar behavior that has been seen in laddcrsi^ 
and square antiferromagnets^ with DM interactions. 

We also studied the effect of the anisotropics on the 
excitation spectra. For that purpose, we have used the 
"bond-wave" theory, that is based on the bosons rep- 
resenting the entangled states of the dimers. In zero 
field we have recovered the momentum-dependent split- 
ting of the triplet states, in accordance with the neutron- 
scattering experiments. Furthermore, we have also re- 
covered the experimentally measured ESR spectra for a 
physically reasonable set of parameters. In this respect, 
we note the foUowings when comparing to the usual per- 
turbational approach, where a single triplet excitation 
propagates^: (i) in order to describe the spectra in the 
field that is parallel to the CuBOs-plane the inclusion 
of all the triplet states is needed; (ii) the finite-field gap 
in the perturbational approach is proportional to the D, 
while in the bond wave approach it is proportional to 
•\/D, thus a smaller value of D can already lead to ob- 
servable effects. 



Regarding the ESR line intensities, we have found 
that in the high-symmetry case, the weight in the spin- 
structure factor appears in the 0(2) phases and is con- 
centrated on the line that is split from the S*^ = triplet 
excitation, and, loosely speaking, follows the "para- 
magnetic" ESR line. In the low-symmetry case, the 
anisotropics make the dimer singlet-0(2) quantum criti- 
cal point a crossover, and the weight accordingly appears 
at the energies of the order of the zero- field singlet-triplet 
gap. This large weight is clearly observed in the ESR 
spectra in the 15-25T range at TOOGHz and above. 

Finally, let us mention the main disadvantage of the 
bond-wave method. Namely, keeping only the quadratic 
terms in the boson operators, the dispersion due to J' is 
not taken into account, and similarly, the J' term does 
not decrease the singlet-triplet gap from its bare value 
A = J (for that we need to go to higher orders in the 1/M 
expansion and to introduce terms with four bosons). In 
other words, the situation is in this respect similar to the 
perturbational approaches that start from the decoupled 
dimcrs. 
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Appendix A: bond-wave theory 

1. General case 

While it is a wclknown procedure, for completeness we 
present how to diagonalize a quadratic form of bosonic 
operators of the form 

where B = B'^ is & hermitian dxd matrix and A = is 
a symmetric d x d matrix, so that the whole Hamiltonian 
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K. is hermitian (however, it is not normal ordered). The 
a^^ denotes a vector of d bosons (a|, Cj, . . . , a^), similarly 
a = (oi, 02, . . . , ad)- We also assume that the 2d x 2(i 
matrix in Ea. (jAip is positive definite - this ensures that 
all the eigenvalues associated with creation operators are 
positive. 

Applying the dO/dt = i[JC,0\ equation of motion of 
an operator O to the bosonic operators a^^ and a we get: 



(A2) 



Our aim is to find a suitable linear combination of a and 
a^^ operators that are energy eigenstates of the JC. This 
is achieved by solving the 





(A3) 



eigenvalue equation. The equation above has a particular 
property: for each (uj,Vj) eigenvector with eigenvalue 
LOj > the (v* , u* ) is also an eigenvector with eigenvalue 
—LOj. We associate the eigenvectors with positive eigen- 
values with creation, and with negative eigenvalues with 
the corresponding annihilation operator: 



and C = ensures that the whole Hamiltonian ma- 
trix is hermitian, moreover we assume it to be positive 
semidefinit. 

Using the equation of motion technique, we get the 
equivalent of the Eq. (|A3p for this case: 



(A7) 




We define the creation operators and energies as 
Oj = Uj • a^ + Vj • b and ujj = w' 



when ujj > 0, and 



a 



^ V* • b'f + u* • a and uj. 



(A8) 



(A9) 



if Wj < 0, so that [ajjaj,] = djj' the commutation rela- 
tions are satisfied. Eventually we arrive at the 

IC' = Y, + ^) + ^ (Tr B - Tr C) (AlO) 

diagonal form. 



a, 



(A4a) 
(A4b) 



Appendix B: Keeping |s) and only 



SO that the [aj,a||,] 



~ "3,3 

filled. With this choice [/C,a| 



-ojjaj holds and /C takes the 



commutation relation is ful- 

Wj-aj and [/C,aj] = 



diagonal form. 



a]a, + - 



(A5) 



In the case of the field parallel to the z axis it is a usual 
practice to keep only the low lying singlet and 5*^ = 1 
triplet (the component aligned with the field) state of a 
bond. Here we are interested in the behavior of the gap 
close to the critical field. For that reason, we restrict the 
discussion to the dimcr-singlet and the 0(2) [6*4] phase 
for D = and the Zi [I?2d] phase in finite D case. 

As a first step, we define the following rotated boson 
operators 



2. Restricted case 

In certain cases the quadratic term is simpler and can 
be written as 



/C' 




(A6) 



where the a^^ = (a|, ai^, . . . , a\,) and b = (61, 62, • ■ • , &d")j 
with d = d' + d" , and the a and h bosons commute among 
each other. This happens frequently for the bosons in the 
momentum representation, as the momentum conserva- 



tion allows only terms of the type oLa 



and aj^a_j. (i.e. 



terms like a\M._-^^ and a-^a-^ are miss- 



ing), and we can associate the Ok and a_\^ bosons with 
a and h bosons in Eq. (|A6|) . respectively. The B — 





= cos 


f^i(k) + 


sin|e^nt^(k). 


(Bla) 


4(k) 


= sin 


f 4(k) - 




(Bib) 


4(k) 


= cos 


f4(k)- 


^sin|e*n^,B(k), 


(Blc) 


tl(k) 


~ sin 


f4(k)- 


icos|e"^4 B(k) , 


(Bid) 



so that the variational wave-function that comprises the 
above mentioned phases is given by 1^')^ ~ |0) ^-nd 

I^)b ^ I*-*)' ^'^'^ phase iy9 = for convenience 

[see also Eq. ((38|) for comparison] . The expectation value 
of the Hamiltonian is then given by 



En — 



— J ( - + COS a 



J' 

— (1 - cosa)^ 
2 ^ ^ 



D'j^ sin^ a 



■—r=D sina — gzhz^l — cos a) 
v2 



(B2) 
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Here we introduce the D = D — 2gshz, as in this section 
the D and the appear in this combination only. Mini- 
mization procedure involves solving a quartic polynomial 
equation that is tedious. Instead, we concentrate on the 
case when the anisotropy terms D is small. 

We also need the bond-wave Hamiltonian. For that 
wc introduce 



symmetric and antisymmetric combination of the rotated 
triplet operators that reduce the size of the matrices in 
the Hamiltonian. Expanding in powers of A/, we get H = 
EoM'^+M'H2 + - ■ ■ , where we omitted higher order terms 
in 1/M. The bond wave Hamihonian H^^) ^ -^(2) +7^(2) 
is given as 



4(k) = i=[4(k)±4 



(B3) 



a±{b + 2D'^)ji 
671 



f±(k) 



a±(6 + 2D^)7i / \ 4(-k) 



(B4) 



where the a can be conveniently expressed as 

J 



and 



= b-Eo-h,-- 



6= ^(J'-2D^)sin2a 



(B5) 



(B6) 



while 7i is defined in Eq. ([5^ . The Bogoliubov transfor- 
mation yields the 



L0± = J(a±2i^^7i)(a±267i±2i?^7i) . (B7) 



1. High symmetry case 



The minimal energy for Z) = of Eq. (|B2[) is achieved 



for 



1 h;, < hci , 

cosao(2) ^ <^ "'-'h.'^XT'' h,,<h,<h,2, (B8) 



-1 



This solutions correspond the the dimer-singlet, 0(2)[5'4] 
and the fully polarized phase, respectively. Note that 
the one-half magnetization plateau is missing - the form 
of the chosen wave function does not allow for the Z2 
breaking. 



Ki = J-2\D'^^\ , 

K2 = J + 2J' + 2|Dl| . 

The variational energy of the unit cell is then 



En — 



_ 3J 

3J _ (h,-h^if 
2 h^2-hci 



i + 2J' - 2g,h, h,2 < K 



(B9a) 
(B9b) 



(BIO) 



-3 




h^/J 



FIG. 15. (a) The lowest lying branch of the excitation spec- 
trum at the k = is shown when keeping 4 (thick) and 2 
bosons (thin lines) per dimer for J' = 0.6 J, = O.IJ and 
different values of D. (b) The bond wave spectrum has a dip 
at the hci = 0.8 J critical field. The dotted line is the approx- 
imation from Ref. [S^, the dashed line and the circle are the 
approximations given by Eqs. (|B17P and (|B19|) . respectively. 



of the 0{2) phase compared to the case when we keep all 
the four state of a dimer [see Eq. and the bound- 

aries overlap only in the limit of small \D'_^\ values, when 
the critical field is close to J. 



It turns out that the boundary between the dimer sin- 
glet phase and the 0{2) phase is shifted to the expense 



Now, let us turn to the excitation spectrum. In the 
dimer-singlet phase a = J — and 6 = in Eq. (|B4p , 
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(2) 

SO that the H^j. matrices arc actually diagonal, 

,(2) 



^r = E^±(k)4(k)f±(k), (Bii) 



with the excitation energies 

^± (k) = J — hz ± 2D'j^ cos ^ cos ^ 



(B12) 



This is the same as the small D'j^ / J limit of the disper- 
sions given by Eqs. ([55l) and ([56|) . when we kept all the 
four bosons per dimer. 

We can also write the bond-wave Hamiltonian in the 



In the case of small D we include the first order correction 
in D to the a given by Eq. (|B8p . 



cos a = \ 1 



D 



yj2(hz - hci){hc2 - hz) ^ 
and we end up with 



5 "0(2) 



^1/2 



{hz - hci)(hc2 - hz) 



1/4 



(B16) 



(B17) 



= E 



J-hz 2Z?l7i\ /<^(k) 



t^k)/ \2Z?l7i J-hz \tB(k) 



(B13) 

form: here we recognize, up to phase factors, the upper 
left corner of the 4 by 4 matrix in Eq. ((5T|) . 

In the 0(2)[S'4] phase the a = J — hd — —2D'j_ and 



^ _ {hz - hci){hc2 - /tz) 
and from Eq. (jB7p wc get 



(B14) 



= (k) = 2Wl Tcos-cos- 



x W D^L ( ^1 =F (^1 + ^) cos ^ cos ^ )(B15) 



We note that w^(k) as k 0, thus it becomes the 
Goldstonc mode associated with the continuous symme- 
try breaking in the 0{2) phase. 



2. Low symmetry case 

In the presence of the D anisotropics the oj^ Goldstonc 
mode acquires a finite gap in the presence of anisotropics. 



in the leading order in D. This approximation is shown 
with dotted line in Fig. [TST b). It clearly fails as ft- — > h^, 
as in the limit ao(2) — ^ the Eq. (jB16P is not valid 
any more. Instead, at the critical field hd and in the 
a — i' limit the energy expression Eq. (|B2p simplifies 
considerably, in leading order 



a — 



and for the gap we get 



2l/6£,l/3 

(J' - 2D'^)^/3 ' 



73^2/3 (J/ „ 2D\ 

22/3 



^l/3 



(B18) 



(B19) 



This approximation is shown with circle in Fig. llSf b). 
Wc find that on the boundary between the dinicr-singlet 
and the 0(2) phase the gap closes faster than D^^^, 
namely with a power 2/3. Such a behavior at the quan- 
tum critical point has been discussed for quantum anti- 
ferromagnets in Refs. [i^ and [i^. We also note that 



the perturbational y {h — hcif + D"^ result of Ref. [38[ 
docs not capture the quantum fluctuation effects close to 
the critical field hd- 
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